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' Abstract 

The present study investigates the interaction of an equidistant three-level atom 
ly-j | and a single-mode cavity field that has been initially prepared in a generalized 

coherent state. The atom-field interaction is considered to be, in general, intensity- 
dependent. We suppose that the nonlinearity of the initial generalized coherent 
state of the field and the intensity-dependent coupling between atom and field are 
distinctly chosen. Interestingly, an exact analytical solution for the time evolution 
of the state of atom-field system can be found in this general regime in terms of 
the nonlinearity functions. Finally, the presented formalism has been applied to 
a few known physical systems such as Gilmore-Perelomov and Barut-Girardello 
coherent states of SU(1, 1) group, as well as a few special cases of interest. Mean 
photon number and atomic population inversion will be calculated, in addition to 
investigating particular non-classicality features such as revivals, sub-Poissonian 
statistics and quadratures squeezing of the obtained states of the entire system. 
Also, our results will be compared with some of the earlier works in this particular 
subject. 

Keywords: Atom-field interaction; Nonlinear Jaynes-Cummings model (JCM); Non- 
linear coherent state; Nonclassical state. 

Pacs: 42.50. Ct, 42.50.Dv, 42.50.Ar, 42.50.-p 
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1 Introduction 



Coherent states (CSs) play an important role in quantum optics and modern physics 
[H EJ |3]. Along different kinds of generalization of coherent states [5], nonlinear CSs [6] 
or /-CSs \a, f) [7J have been introduced and attracted much attention in recent decade 
[8]-[T3] with the Fock space representation as [7J: 

/ oo |q,| 2 ™ \ _1 ^ 2 °° n 

la, /) = > ... —, ,,,. > — p= In), a = \a\e tip , (1) 

where [/(^)]1 = f(l)f(2)...f(n). According to this formalism, /-deformed annihilation 
and creation operators, respectively are defined as A = af(n) and so = f^{n)a) where 
a, and n = a^a are bosonic annihilation, creation and number operator, respectively. 
Here, the (real) intensity-dependent function f(n) is responsible for the nonlinearity of 
the states. On the other side, Jaynes-Cummings model (JCM) is the simplest nontrivial 
example of the atom- field interaction, i.e., a two-level atom and a single-mode radiation 
field [14J. As a few recent works in this topics see Refs. [T51 [TB] • The dynamical behaviour 
of the nonlinear atom-field interaction in the presence of classical gravity using the non- 
linear coherent states approach discussed in [17] . Buzek generalized the JCM [18] and the 
atom-field coupling is considered to be intensity-dependent and supposed that the cavity 
field be in the Gilmore-Perelomov (GP) nonlinear CS of SU(1, 1) group [19]. The author 
showed that, the revivals of the radiation squeezing are strictly periodical for any value 
of initial squeezing parameter [IB]. More recently, Koroli et al [20] studied the interaction 
of an equidistant three- level atom (ion), whose dipole moment matrix transition elements 
between the adjacent atomic levels are different, with the GP CS of SU(1, 1) group. They 
showed that, in the three-level model with the intensity-dependent coupling, the exact 
periodicity of the squeezing revivals is violated. 

In the present paper, due to the great interest in the atom-field interactions in the 
quantum optics, we regard the same configuration of three-level atom (ion) has been 
considered in [20] , however we generalize the initial state of the field to the " nonlinear CS" 
with arbitrary "nonlinearity function" f(n). In our formalism, the atom-field coupling 
is also considered to be intensity-dependent, characterized by a function g(n), which 
is generally different from f(n). While in Refs. [HJ [20] the authors have concerned 
with special nonlinearities, our presentation deals, in principle, with two distinct general 
nonlinearity functions and interestingly, the exact solution is also obtained for the time 
evolved entire states of the system. 

After finding the explicit solution of the state vector, which evidently depends on f(n) 
and g(n), we will treat the time evolution of the mean photon number and atomic popu- 
lation inversion, as well as a few non-classicality features such as sub-Poissonian statistics 
and squeezing of the quadratures of the field. Therefore, in contrast to the approach of 
[20] , that investigated only a special system, i.e., Holstein-Primakoff SU (1,1) coherent 
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states as filed state and y/n as the coupling between atom and field, our formalism can 
be applied to the intensity-dependent interaction between the same three-level atom with 
arbitrary cavity field that is initially prepared in a generalized CS with known "non- 
linearity function /(n)". Interestingly, apart from the mentioned advantage, as we will 
observe, our presented formalism allows one to work with the general case in which two 
different nonlinearity functions exist, one indicates the initial state of the field (f(n)), 
and the other to the atom- field intensity dependent coupling (g(n)). Although, it can 
simply recover the work of Buzek [18] and Koroli et al [20], as some special cases. In the 
continuation, the formalism has been applied to two classes of generalized CSs, i.e., two 
distinct representations of SU (1,1) CSs. In addition, in opposite of the work has been 
done in [20], in which the authors have taken f(n) = 1/ 'y/n and g(n) = \/n, we choose 
f(n) = g(n) = \/n, i.e., the nonlinearity of the initial state is taken to be the same as the 
intensity-dependent coupling function. Finally, in each of the considered cases, we dis- 
cuss the numerical results have been shown in several figures in detail and compare with 
earlier works. Along this procedure, some new physical features reveals, which may be 
emphasized. For instance, a typical collapse and revival phenomenon (as a pure quantum 
mechanical feature) in physical quantities such as mean photon number, atomic popula- 
tion inversion and Mandel parameter, and the appearance of some non-classicality signs, 
are specific aspects may be highlighted. 



2 Hamiltonian of atom-field system 

We consider the interaction of an equidistant three-level atom in a cascade configuration 
with different dipole moment matrix transition elements between the adjacent atomic 
energy levels (shown in Figure 1) with a quantized single-mode cavity field of frequency 
co. The states \g), fa) and fa) are respectively denoted the ground, first and second 
excited states. Only the atomic transitions between \g), fa) and also fa), fa) are dipole 
allowed, but not between \g), fa)- The Hamiltonian of such an atom-field system is given 
by: 

H = H + H U (2) 

where H is the sum of the Hamiltonians of atom and field, i.e., H = huoS z + hua'a, 
S z — fa)fa\ — \g)(9\i denotes the atomic population inversion operator. The interaction 
Hamiltonian between atom and field in (j2j) in the dipole and rotating wave approximation 
is given by: 

E x = ^Ai(a t |e 1 )(e 2 | + a|e 2 )(e 1 |) + h\ 2 (a ] \g) (e x | + afa)(g\), (3) 

where Ai,A2 are the atom-field coupling constants. A useful approach to the atom-field 
interaction problem may be found in the "interaction picture" [21J. The Hamiltonian 
describing the interaction between the mentioned atom-field, in the interaction picture, 
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is described by 

Hj = ^A 1 (at|e 1 )(e 2 |e- iA * + a|e 2 ) ( ei |e iA *) + HX 2 (a^\g)(e 1 \e- iAt + a\e 1 ){g\e iM ) (4) 

where, A = u — oj is the detuning. Now, we suppose that the atom-field coupling is 
intensity dependent, expressed by g(ri), so the intensity-dependent Hamiltonian in the 
interaction picture is given by: 

Hi = nA 1 (i? t |e 1 )(e 2 |e- lA ' + J R|e 2 )(e 1 |e iA *) + hX^R^g)^ 1 ^ + R\e x ) (g\e lAt ) (5) 

where R = ag(n), fit = g(n)a) and g(n) describes the intensity-dependent coupling 
between atom and field. Notice that, R and R) in ([3]) have the same structure and also 
meaning of A and were described at the beginning of this section, except that a different 
nonlinearity function is considered. Our different notation is only for the distinction 
between the nonlinearity of the initial state of the field and the intensity dependent 
coupling function. We have assumed that g(n) is a real well-defined function with no 
singularity. The Hamiltonian in (jHJ) plays a crucial role in determining the subsequent 
dynamics of the quantum states of a variety of (three-level) atom-field systems, some of 
them will be considered in the continuation of the paper. 



3 Atom-field state vector 

Let the atom be initially in the first excited state |ei) and the cavity field is prepared in 
a generalized CS with a nonlinear CS as given in ([1]). The state of the atom- field system 
at t = can be expressed as: 

oo 

= 0)) = | ei ) ® \a,f) = J2 C nh,n), (6) 

n=0 

where C n is determined as the expansion coefficients of the states in (pQ), according to the 
initial CS has been chosen. A useful approach to the atom-field interaction problem may 
be found in the interaction picture |21j. The Schrodinger representation of a state vector 
|\l/(t)), in terms of its interaction picture representation | *&/-(£)), is given by |^(t)) = 
U Q (t)|tfj(t)), where U (t) = exp(^f^). Also, can be obtained from zA^|*j(t)) = 

Hi\^fi(t)). Consequently, by using the above relations, one can straightforwardly find 
the time evolution of the state vector of the coupled atom-field system in the resonance 
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condition A = as follows: 



!*(*)> 



^ y e -iuj {S z +n)t q 



n=0 



x {cos(a/ ng 2 (n) + (3 2 (n + l)g 2 (n + 1) r) \ei, n) 
(3Vn + lg(n + l) 



^ng 2 (n) + f3 2 (n + l)^ 2 (n + 1) 
x sin(v / ns ,2 (™) + /3 2 (n + l)g 2 (n + 1) r) |gr,n + 1)} 



^ ^ e -iwo(S s +n)l q _ 



ng{n) 



^ng 2 (n)+f3 2 (n+l)g 2 (n + l) 
x sm(^/ng 2 (n) + (3 2 (n + l)g 2 {n + 1) r) |e2, — 1), 



(7) 



where /3 = A2/A1 and r = Ait. We will call r which is the scaled time briefly as time. 
It is worth noticing that n which is appeared in the exponential functions in (JJJ) has the 
operational role (number operator), the same as the role of S z . Anyway, this solution 
is very general and gives the time evolution of a system that is involved a single-mode 
cavity field, has been prepared in a generalized CS, coupled with an equidistant three- 
level atom in a cascade configuration with different dipole moment matrix transition 
elements between the adjacent levels (Ai 7^ A2). Two specific cases can be recovered: 
(1) when A2 — > (/3 — > 0), it is equivalent to the single two-level atom; and (2) when 
A2 — > Ai ((3 — > 1) it corresponds to an equal dipole moment matrix transition elements 
between the adjacent levels. This case is indeed equivalent to a pair of indistinguishable 
two-level atoms [22J. The states of the atomic pair can be described in the three-level 
states representation: (i) the ground state is equivalent to the case in which both atoms 
of the pair are in the ground state; (ii) the first exited state describes the case in which 
one atom of the pair is in the ground state and another atom is in the exited state, and 
(iii) in the second exited state, both atoms of the pair are in the exited state. In addition 
to the generality of our presented formalism, another advantage of our formalism is that 
it contains two distinct nonlinearity functions, i.e., f(n) corresponding to initial state of 
the field and g(n) which determines the intensity-dependent coupling between atom and 
field. So, we can easily recover the results in [20], if we take f(n) = l/y/n and g(n) = \/n 
and also in [18] . if we take f(n) = l/y/n, g(n) = \fn with /3 = 0, as some special cases. 
Also, we can work with other different possibilities, especially the case g(n) = f(n) = 1 
and g{n) = f{n) ^ 1, investigate the output results and compare with [20] . 
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4 Quantum statistics and non-classicality of the in- 
troduced state 

After determining the state vector of the atom-field system at any arbitrary time obtained 
in (JTj), we are able to investigate mean photon number, atomic population inversion, 
Mandel parameter and squeezing parameters as some non-classicality criteria. 



4.1 Mean photon number 

Using the explicit form of atom-field system in (j7|), the mean photon number can 

be readily found for any system, with arbitrary f(n) and g(n) as follows: 

CO 

( n ) = \C n \ 2 n cos 2 ( y/ng 2 (n) + (3 2 {n + l)g 2 (n + 1) r) 



n=0 

oo 



+ l^l 2 sm 2 {y/ng 2 (n) + (3 2 {n + l)g 2 {n + 1) r) 

n=0 

n(n - l)g 2 (n) + (3 2 (n + l) 2 g 2 (n + 1) 
X ng 2 (n) + (3 2 (n + l)g 2 (n + 1) ' ^ } 

4.2 Atomic population inversion 

Atomic population inversion, as the expectation value of S z , is defined as 

(s z ) = m)\s z \nt)) = \(e2\nt))\ 2 - mm)\ 2 , (9) 

where |(e2|^/(t))| 2 and |((7|\l/(t))| 2 are the probabilities of the presence of the atom in \e2) 
and \g) states, respectively. With the help of in (JTj), (S z ) can be written as: 

00 

(S z ) = \ C n\ 2 sm 2 (^/ng 2 (n) + (3 2 (n + l)g 2 (n + 1) r) 

n=0 

ng 2 (n)-/3 2 (n + l)g 2 (n + l) 

ng 2 (n) + f3 2 (n+l)g 2 (n+lY 1 } 

As a special case, in the absence of the radiation field inside the cavity, namely, the cavity 
field being in the vacuum state (n = 0), the quantity (S z ) is given by: 

(^) n=0 = -|C | 2 sin 2 (/3^(l)r). (11) 

Recall that, the atom is initially in \e%). Therefore, in the absence of a driving field, the 
atom in the lower state |ei), cannot excite to the upper state |e 2 ), so (|(e2|^ f (t))| 2 ) n=0 = 
and from one has 

(^)„=o = -(l(5l^)>| 2 ) n= o- (12) 
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By the comparison of ( ITT]) and (fT2|) . one can arrive at (| (g\^/(t)) | 2 ) n=0 = |C*o| 2 sin 2 (/3 r), 
which means that in the fully quantum theory of radiation, transition from the upper state 
to the lower state in the vacuum of the field becomes possible, known as the spontaneous 
emission. It is to be noted that this result cannot be predicted by semiclassical radiation 
theory. 

4.3 Mandel's Q-parameter 

To examine the statistics of the states, Mandel's Q-parameter is widely used, characterizes 
the quantum statistics of the states inside the cavity. This parameter has been defined as 
Q = (w2) ( - (n)2 - 1 [23], where (n) obtained in ® and (n 2 ) may be calculated as follows: 

oo 

(n 2 ) = J2 \°n\ 2 {n 2 cos 2 (vWW + P 2 {n + l)g 2 {n + 1) r) 

n=0 

n(n - l) 2 g 2 {n) + (5 2 {n + l) 3 g 2 (n + 1) 
+ ng 2 (n) + (3 2 (n + l)g 2 (n + 1) 

x sin 

(vWH + P 2 (n + l)g 2 (n + 1) r)\ . (13) 

The Mandel's Q-parameter obviously depends on the particular choice of f(n) and g(n). 
It is well-known that Q is positive for classical states (super- Poissonian), negative for 
non-classical states (sub- Poissonian) and vanishes for canonical CSs (Poissonian). 

4.4 Squeezing parameters 

Also, we will investigate the squeezing properties of the quadratures of the field. For this 
purpose, we introduce field quadratures as Xi = , X 2 = A ~ A , where A and are 
the operators A = a e luJt ,A^ = e~ luJt . To study the squeezing properties, we introduce 
the squeezing parameters: 

^■(r)=4((AX,) 2 )-l, (14) 

where Sj(r) corresponds to squeezing effect in Xj and satisfies the inequalities —1 < 
Sj(r) < 0. Obviously, to preserve the Heisenberg uncertainty relation, when Si^) is 
negative, S^Si) should be positive. To calculate the parameters Sj(r) numerically, one 
has to find the mean values of the operators: a,a 2 ,a^ and a) 2 . The following results are 
easily obtained: 

(a) = e-^-^B^r), (a 2 ) = e" 2 ^'-^ B 2 (t) , (15) 
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where 



oo 



?1=0 



x cos(^/ng 2 (n) + f3 2 (n + l)g 2 (n + 1) r) 
x cos(a/(t2 + l)g 2 (n + 1) + (3 2 (n + 2)g 2 (n + 2) r) 
sin(A/ng 2 (n) + /3 2 (n + l)g 2 (n + 1) r) 



+ 



^ng 2 (n) + /3 2 (n + l)g 2 {n + 1) 



sin 
x — 



x 



sin(v/(n + l)# 2 (n + 1) + (3 2 {n + 2)# 2 (n + 2) r) 
V^ra + l)# 2 (n + 1) + f3 2 (n + 2)# 2 (n + 2) 

x v^9(n + l)N(n) + ^ 2 (n + %(n + 2)]} (16) 



OO 

B 2 (t) = e- 2 ^J2 C nCn + 2 
n=0 

x cos(^/ng 2 (n) + (3 2 (n + \)g 2 (n + 1) r) 
x cosfV (n + 2)# 2 (n + 2) + /3 2 (n + 3)# 2 (n + 3) r) 
sin(A/ng 2 (n) + /3 2 (n + l)g 2 (n + 1) r) 
^ng 2 (n) + /3 2 (n + l)g 2 (n + 1) 
(y/(n + 2)g 2 (n + 2) + f3 2 (n + 3)g 2 (n + 3) r) 
v/(n + 2)# 2 (n + 2) + /3 2 (n + 3)# 2 (n + 3) 
\/n + l\/n + 2[ng(n)g(n + 2) 
+ /3 2 (n + 3Mn+l)s(n + 3)]}. (17) 

Clearly, (a*) = (a)t and (a 1 " 2 ) = (a 2 )*. It should be noticed that Bi(t) and B 2 (r) are real 
values. Finally, using the above expressions we arrive at 

S 1 (t) = 2[B (r) - B 2 {r)\ + 4 cos 2 (y?) [B 2 (t) - B\{t)\, (18) 

S 2 (r) = 2[B (t) - B 2 {r)\ + 4 sin 2 (y?) [B 2 (r) - B 2 (r)}, (19) 
where we have set B (t) = (n). 

5 Some physical realizations of the formalism 

In this section, we firstly consider the special simple case f(n) = g(n) = 1 and then 
nonlinearity functions of CSs of 577(1, 1) group are considered. 
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5.1 The special case: f(n) = g(ri) = 1 

In this subsection it will be constructive to consider canonical coherent state as the initial 
state of the field and suppose that the coupling between (the three-level) atom and field 
be independent of intensity, namely f(n) = g(n) = 1. The quantities that have been 
presented in section 4, are plotted versus time for this special case. From figures 2-4 it is 
seen that mean photon number, atomic population inversion and Mandel parameter have 
nearly complete collapses and revivals, at least at intermediate times. This quantities 
collapses to (n) ~ 63.5, (S z ) ~ 0.5 and Q ~ 0.01, respectively. Atomic population 
inversion is positive. This means that the probability of the presence of the atom in |e2) 
is greater than being in \g). Mandel parameter is negative in some time intervals, namely 
the state of the system possesses sub-Poissonian statistics. Squeezing parameters versus 
time is shown in Fig. 5. There is no squeezing in X\ and nearly in X2, too if one ignores 
the very weak squeezing that may be seen about r = 45. 

5.2 Nonlinear CSs: CSs of 577(1, 1) group 

At this stage of the paper, as an example of initial nonlinear CS, we will consider the 
Holestein-Primakoff single-mode realization of SU(1, 1) Lie algebra. Before we proceed, 
it is worth noticing that, the equivalence of the discrete series representation of SU(1, 1) 
state space {\k, n)}^L , with k — ~, 1, |,2, ... , and the harmonic oscillator Hilbert space 
{\n)}^ = Q is illustrated in [21]. Based on this recognition, SU (1,1) CSs have been well 
established as nonlinear CSs by Ali et al [25J. 

We briefly introduce the GP CSs of SU(1,1) group (sometimes have been called 
Klauder-Perelomov CSs [21 E]). These states |20 are defined in the interior of the unit 
disk in the complex plane, centered at the origin. The nonlinearity function corresponding 
to these states is deduced as fcp(n, k) = l/y/n + 2k — 1 [U [HI [27]. In this subsection, 
we first take f(n) = g(n) = fcp(n,K) in our further numerical calculations. Different 
quantities mentioned in section 4 associated to the atom- field state (J7J), have been plot- 
ted versus time in figures 6-9 with fixed parameters k = 3/2, \a\ = 0.9 and = 0.01. 
Figure 6 deals with the mean photon number. A typical fractional collapses and revivals 
are visible from the figure 6a. This figure indicates that the envelope of the oscillations 
fractionally collapses to a fixed value ~ 12.3 and as time goes on, the collapsed mean 
photon number is partly revived. The maximum amplitude of the oscillations occurred 
at t = decreases with time, and the duration of the oscillations varies, irregularly. Our 
aim for showing the figure in a short time interval is to explain the details of the varia- 
tion of mean photon number in a more apparent fashion, particularly in relation to next 
figure. Atomic population inversion is shown in figure 7. The positivity of this quantity 
at all time means that the probability of the presence of the atom in the state \e^) is 
larger than the probability of being in the state \g). Similar to mean photon number, 
fractional collapses and revivals are observable from figure 7a. It is to be noted that, in 
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this case the atomic population inversion collapses to a fixed value ~ 0.5. Figures 6b 
and 7b are the same as figures 6a and 7a, respectively, with the same chosen parameters, 
except that the interval of time are restricted to 20. A comparison between the figures 6b 
and 7b shows that they both obey nearly the same general pattern and makes one sure 
that the variations of mean photon number and the atomic population inversion are in 
opposite directions, as is expected. Figure 8 displays Mandel parameter as a function of 
time. The supper- Poissonian behaviour is observed from the figure, together with a (typ- 
ical) collapses and revivals. This phenomenon occurs for the system under consideration, 
although not in a clear and regular manner. Indeed, Mandel parameter oscillates such 
drastically with time that the quasi-chaotic behaviour is revealed. In figure 9 the curves 
of squeezing parameters Si and S2 are shown for tp = ir/2. Only the negativity of Si, 
which indicates squeezing effect in X\ quadrature, is revealed in a finite interval of time. 

At this stage, in contrast to Buzek [18J and Koroli et al [20J that took f(n) = 
fcp{n, 1/2) = l/y/n and g(n) = \/n, we will set fscin) = g(n) = y/n and follow our 
numerical calculations. The nonlinearity function \fn can be associated with a particular 
case of the BG state when one chooses k = |. These states are defined in the whole 
of the complex plane. We emphasize that as in the GP CSs, we consider the Holestein- 
Primakoff realization of SU(1, 1) Lie algebra of BG states. Barut-Girardello (BG) CSs 
of SU(1, 1) group [28] are established as the dual pair of GP CSs [25]. The nonlinearity 
function of these states is obtained as f B c(n, k) = \/n + 2k — 1 [TT], [25J [2?]. It is easy to 
check that the operators K_ B q = afBa{n), K+,bg = f BG {n)a\ K^ BG = ±[K_ BG , K +jB g] 
satisfy the commutation relations [K- )BG , K +)BG ] = 2K 0tBG , [K Q ^ BG , ^±,bg] = ±K± >BG . 
where K = n + k. 

Anyway, choosing f{n) = g{n) = f B G(n, 1/2) = y/n, our calculated results have been 
displayed in figures 10-13. The mean photon number and atomic population inversion 
versus time for the corresponding atom-field states have been shown in figures 10 and 
11, respectively. In the range which the atomic population inversion gets negative values, 
the probability of the presence of the atom in \g) is indeed more than |e2). Observing 
figures 10 and 11, it will be clear that, any increase (decrease) in population inversion 
is simultaneous with a decrease (increase) in mean photon number. Also, there are two 
distinct oscillatory behaviour in these two quantities, one (small) is within the other 
(large). This feature was not appeared in [20], when they took /(n) 7^ g(n) (figures 1 
and 2 of Ref. [20]). While there are some jumps in the mean photon number (which 
simultaneously accompanied by a sudden decrease in the atomic population inversion), 
according to our results small jumps and decrease within great jumps and decrease will 
be revealed. It is worth mentioning that, we continued our numerical results (have not 
shown here) with calculating the mean photon number and atomic population inversion 
in a wide interval of time and observed that, for both values of the considered 0, the 
quantities are nearly (not exactly) periodic. But, the period of time for the case (5 = 0.01 
is larger than the case = 0.1. In figure 12, the Mandel parameter has been displayed as 
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a function of time. The results indicate that for (3 = 0.1 this parameter is almost periodic, 
while its maximum amplitude decreases with time. The sub-Poissonian behaviour also 
occurs in wide intervals of time. For large times, full sub-Poissonian statistics will be 
revealed for such system, in both of the presented cases. It is worth noticing that, this 
non-classicality sign has not been observed in |20J. In figure 13, the curves of squeezing 
parameters repeated regularly, but not exactly. Only for small intervals of time, squeezing 
effect is occurred in X\ or X 2 quadrature. 

Also a specific case can be regarded, i.e., when f(n) = 1, which means that the 
initial field is in the standard CS and the interaction between atom and field would be 
intensity-dependent. For instance, we consider g{n) = fcp(n,K) = l/y/n + 2k — 1. Our 
numerical results have been displayed in figures 14-17. In figure 14 the mean photon 
number for the atom-field state ((?]) associated with these particular functions is shown. 
Atomic population inversion is plotted versus time in figure 15. A careful observation on 
the figures 14, 15 leads one to conclude that, when the mean photon number of the field is 
increased (decreased) due to emission (absorbtion), the probability of the presence of the 
atom to be in the \g) (|e2)) state increases, consequently the atomic population inversion 
is reduced (increased). Figure 16 shows the Mandel parameter, which implies that it is 
negative in a wide range of time, so that sub-Poissonian behaviour occurs, repeatedly. 
As shown in figure 17, in small regions of time the squeezing effect is observed weakly in 
X\ or X 2 quadratures, separately. Obviously, when SxiS-i) is negative, Si{S\) have to be 
positive. 

Summing up the above presented results, two remarkable points may be offered. 
Firstly, the represented numerical results plotted in figures 2-5 for the case f(n) = g(n) = 
1 show that the collapses and revivals occur nearly regularly (but not exactly) relative to 
the next ones which contain some kind of nonlinearities. Indeed, in the latter cases the 
chaotic behaviour will be revealed due to the presence of the nonlinearities. In general, 
the Mandel parameter and squeezing effect have not an exact regular periodicity, specif- 
ically relative to two-level atoms. In addition, as it is expected, the variations of mean 
photon number and the atomic population inversion are in opposite directions. As we 
have explained, this is consistence with the physics of the considered interaction. Also, 
the fractional collapses and revivals phenomenon, as a well-known non-classicality sign, is 
seen in the mean photon number, atomic population inversion and Mandel parameter in 
the two groups of figures 6-8. The negativity of Mandel parameter in a wide range of time 
is observed in figures 12 and 16. The latter effect, will become more important, if one 
recalls that the squeezing signature of the field is revealed only in a small finite interval 
of time (see figures 9, 13, 17). 
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6 Summary and concluding remarks 



In summary, we have considered the interaction of an equidistant three-level atom in a 
cascade configuration with different dipole moment matrix transition elements between 
the adjacent atomic levels, and a single-mode cavity field that is initially prepared in 
a generalized CS. It should be emphasized that, our formalism is presented in a very 
general regime, because it contains two nonlinearity functions: f(n) which characterizes 
the initial state of the field and g(n) which determines the intensity-dependent coupling 
between atom and field. Particularly, f(n) = g(n) = 1 is equivalent to the case in which 
the initial state of the field is standard CS and the interaction between atom and field 
be independent of intensity. Interestingly, we have presented a closed analytical solution 
for such a non-trivial problem. Then, as some physical appearances of the proposed 
structure, we have investigated the mean photon number, atomic population inversion, 
Mandel parameter and squeezing parameters for GP and BG CSs of SU(1, 1) group as 
well as some special cases. A few points are remarkable, regarding the presented results. 

• Unlike the reported work by Buzek in [18] which considered the two-level atom, we 
have not observed the exact regular periodicity of the squeezing parameters, in nei- 
ther of the chosen nonlinearity functions. We also calculated the Mandel parameter 
and find that the exact periodicity which occurs in two-level atom violates. These 
phenomenon are consistence with the reported results in [20]. We investigated and 
examined this observations for various cases, either with the same nonlinearities 
(f(n) = g{ri)) or with different nonlinearities (f(n) ^ g{n)). 

• The variation of mean photon number and atomic population inversion are in op- 
posite directions, which is an expected result, in view of the physics of the sys- 
tems under the considered interaction. Also, for the special case fcp( n ) = g{n) = 
1 / \/n + 2k — 1 in the second group of figures, the fractional collapses and revivals 
phenomenon in the mean photon number and the atomic population inversion, as 
well as the Mandel parameter are new features of our proposal, may be highlighted. 

• Our results confirm that only in the fully quantum theory of radiation, the sponta- 
neous emission (transition from the upper state to the lower one) in the vacuum of 
the field becomes possible. 

• Comparing the figures 2-5, represented the numerical results have been plotted for 
the case f(n) = g(n) = 1, with next figures one can see that the collapses and 
revivals occur more regularly for the former ones relative to others which contain 
some kind of nonlinearities, either in nonlinearity function of initial state or the 
intensity-dependent coupling. Indeed, in the latter cases the chaotic behaviour will 
be revealed due to the presence of the nonlinearities. 
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• Apart from these, the generality and at the same time the simplicity of our proposal 
allows one to apply it to other physical systems with known nonlinearity functions, 
for instance, center of mass motion of trapped ion [6], photon-added CSs [10], q- 
CSs [7], deformed photon-added nonlinear coherent states [29] have been recently 
introduced by one of us and so on. These are straightforward tasks may be done 
elsewhere. On the other side, Roknizadeh et al introduced a Hamiltonian associated 
with a nonlinear oscillator system, based on action identity requirement of nonlinear 
CSs as follows: H = A^A = nf 2 {n) [llj, upon which one obtains the eigenvalues 
e n = nf 2 {n) or equivalently f(n) = • So, obviously the presented approach can 
be easily applied to such one- dimensional solvable systems, too. 
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FIGURE CAPTIONS: 



Figure 1 Schematic diagram of the cosidered three-level atom. 

Figure 2 The time variation of (n) with |a| = 8 and f3 = 0.01 for f(n) = g(n) = 1. 

Figure 3 The time variation of (S z ) with \a\ = 8 and (3 = 0.01 for f(n) = g(n) = 1. 

Figure 4 The time variation of Q with \a\ = 8 and (3 = 0.01 for f(n) = g(n) = 1. 

Figure 5 The time variation of S\ (the solid curve) and S2 (the dashed curve), with 
|ot| = 8, p = 0.01 and <p = vr/2 for f(n) = g(n) = 1. 

Figure 6 The time variation of (n) for different intervals of r, with k = 3/2 with 
|a| = 0.9 and (3 = 0.01 for f(n) GP = g G p{n) = l/Vn + 2/c- 1. 

Figure 7 The time variation of (S z ) for different intervals of r, with k — 3/2 with 
|a| = 0.9 and (3 = 0.01 for / GP (n) = g G p(n) = l/Vn + 2/c- 1. 

Figure 8 The time variation of Q for k = 3/2 with |a| = 0.9 and (3 = 0.01 for 

/gp(w) = 9Gp{n) = l/\ / n + 2n- 1. 

Figure 9 The time variation of Si (the solid curve) and S% (the dashed curve), 
for k = 3/2 with |a| = 0.9, (3 = 0.01 and ip = it/2 for fap{n) = Qcpij 1 ) — 
l/Vn + 2/c- 1. 

Figure 10 The time variation of (n) for /3 = 0.1 (the solid curve) and [3 = 0.01 
(the dashed curve) with |a| = 0.6 for f(n) = g(n) = \fri. 

Figure 11 The time variation of (S z ) for /3 — 0.1 (the solid curve) and (3 = 0.01 
(the dashed curve) with \a\ = 0.6 for f(n) = g(n) = yjn. 

Figure 12 The time variation of Q for (3 — 0.1 (the solid curve) and (3 = 0.01 (the 
dashed curve) with |a| = 0.6 for f(n) = g(n) = y/n. 

Figure 13 The time variation of Si (the solid curve) and S2 (the dashed curve), 
with I ok j = 0.6, (3 = 0.1 and ip = n/2 for f(n) = g(n) = y/n. 

Figure 14 The time variation of (n) for k — 2 (the solid curve) and k — 4 (the 
dashed curve) with |a| = 0.5 and [3 = 0.1 for f(n) = l,gap(n) = 1/Vn + 2k — 1. 

Figure 15 The time variation of (S z ) for k — 2 (the solid curve) and k = 4 (the 
dashed curve) with |a| = 0.5 and (3 = 0.1 for f(n) = l,gap(n) = 1/ \fn + 2k — 1. 

Figure 16 The time variation of Q for k = 2 (the solid curve) and k — 4 (the 
dashed curve) with \a\ = 0.5 and j3 — 0.1 for f(n) = l,gap(n) = 1/ ' \fn + 2/t — 1. 

Figure 17 The time variation of Si (the solid curve) and S2 (the dashed curve), for 
k — 2 with I I = 0.5, (3 = 0.1 and ip = n/2 for f(n) = l,gap(n) = 1/ V n + 2k — 1. 
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